Heavy-Light Fermion Mixtures at Unitarity 
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We investigate fermion pairing in the unitary regime for a mass ratio corresponding to a 6 Li - 
40 K mixture using Quantum Monte Carlo methods. The ground-state energy and the average light 
and heavy particle excitation spectrum for the unpolarized superfluid state are nearly independent 
of the mass ratio. In the majority light system, the polarized superfluid is close to the energy of a 
phase separated mixture of nearly fully polarized normal and unpolarized superfluid. For a majority 
of heavy particles, we hnd an energy minimum for a normal state with a ratio of ~ 3:1 heavy to 
light particles. A slight increase in attraction to kFa w 2.5 yields a ground state energy of nearly 
zero for this ratio. A cold unpolarized system in a harmonic trap at unitarity should phase separate 
into three regions, with a shell of unpolarized superfluid in the middle. 
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Superfluid pairing and the equation of state of cold 
trapped atoms in the unitary regime have recently been 
the subject of intense theoretical and experimental in- 
vestigation [J Q. These systems are closely related to 
strongly interactingfermions in other regimes, such as 
neutron matter [3,13] and dense quark matter [H, and 
hence are useful as prototypes and benchmarks in many 
areas of physics. At unitarity, all physical quantities are 
simply given by dimensionlcss numbers times the rel- 
evant free Fermi gas quantity. Theoretical predictions 
based upon Quantum Monte Carlo (QMC) calculations 
for these dimensionless numbers - including the ground- 
state superfluid energy £ = E sf /Epo w 0.4 [a 0, H, 
the pairing gap r\ = A/Ep as 0.5 @, H, [j|, and the first- 
order phase transition between an unpolarized superfluid 
and a normal state at finite polarization or concentration 

l 1211 - ar e in good agreement 



/n T « 0.44 0, E 



with recent experiments 13, 1J, 1 

An intriguing variation of this problem is pairing be 
tween particles with different masses, which is within ex 
perimental reach (Til . 17 1 and has already sparked con 



siderable theoretical interest 0, [H, H3, HH, [H [H [24, 
25 . [26|. The most promising candidate is a mixture 
of 6 Li and 40 K s-wave Feshbach resonances, for which 
the mass ratio r « 6.5. A heavy-light fermion mixture 
may be more likely to exhibit exotic phases, like Larkin- 
Ovchinnikov-Fulde-Ferrell phases [27j, while for higher 
mass ratios or more attractive interactions Efimov states 
are expected to appear. 

We consider an interaction of the form: 
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where h denotes a heavy particle and I denotes a 
light particle, with a mass ratio r = rrih/mi, and a 
zero-range interaction between light and heavy parti- 



cles with strength tuned to infinite scattering length 
in the unequal-mass pair. Mean-field BCS theory for 
unequal-mass pairing predicts a simple scaling of the 
equation of state in terms of the reduced mass m r = 
mimh/(mi + ran). If we define the average chemical po- 
tential by p, = {nh + mO/2 then fx and the pairing gap 
A remain unchanged in units of the reduced Fermi en- 
crgy E™- = -^-(37r 2 n) 2 / 3 = . where n is the total 
particle density. 

The heavy and light excitation energies naturally de- 
pend upon the masses rrih and mi individually. The en- 
ergies of the heavy and light excitations are: 

m &(i)(fc)-fr(/o(*0 
&h(l)W = 5 h 



+ A 2 (fc), (2) 



where (k) = - Li h[l) . Even so, the average of 

Eh{k) and Ei{k) depends only upon the reduced mass 
m r , as does the gap A(fc). 

There is no a priori reason to believe that the BCS 
results should be accurate. We have performed QMC 
calculations of the homogeneous superfluid phase, ex- 
amining the quasi-particle dispersion as a function of 
the momentum. The methods are those employed pre- 
viously in the equal-mass case [H, Q, using a modified 
Poschl- Teller potential with an effective range of ro/12, 
where 4/37rrg = 1/n. The superfluid and normal phase 
trial wave functions are of the same form as used pre- 
viously, and provide fixed-node upper bounds to the en- 
ergy; the superfluid wave function has been variationally 
re-optimized. Of course, new physics corresponding to 
quite different nodal strucures of more exotic trial func- 
tions (for example LOFF phases) is not excluded. 

For a mass ratio of 6.5, we obtain a ground-state energy 
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FIG. 1: (color online) Quasi-particle excitation energies for 
the two different species. The QMC results are shown as 
diamonds and squares for light and heavy particles, respec- 
tively. The BCS results are shown as lines: the dashed line 
corresponds to light particles, the dot-dashed line to heavy 
particles, while the solid line to the results for the equal-mass 
case. Also shown are the QMC results (circles) from Ref. @] 
for the equal-mass case, as well as the average (triangles) of 
the two QMC sets of points for heavy and light particles. 



£(r = 6.5) = 0.390(5), slightly lower than the £(r = 
1) = 0.41(1) obtained for the same interaction with equal 
masses. The latter extrapolates to £(r = 1) = 0.40(1) at 
zero effective range; we have verified that similar small 
extrapolations are present for unequal masses. The small 
difference (~ 5%) between the ground-state energies from 
r = 1 to r = 6.5 is a measure of the contribution of 
non-zero total momentum pairs in the ground state. In 
QMC calculations the ground state is determined by a 
diffusion algorithm which depends upon the local energy 
and the mass of the particles. For a pure BCS state, the 
local energy H^t/^t is identical over the entire 3N- 
dimcnsional coordinate space for different mass ratios. 
However the total mass of the pairs increases with r, 
resulting in a slightly lower energy as the diffusion of the 
pair centers-of-mass is reduced. 

We have also calculated the quasi-particle excitation 
energies for the light and heavy particles; the results 
are shown in Fig. [TJ The excitation energies shown for 
the light (heavy) particles are calculated by subtracting 
the energy of N/2 heavy and N/2 light particles from 
the energy using a trial function with an additional light 
(heavy) particle in a state of momentum k. Both simula- 
tions are performed for the same volume L 3 = 3ir 2 N/kp. 
To facilitate comparison with the QMC results, the BCS 
lines shown are Eh(k) + fih an d Ei(k) + (ii for the heavy 
and light particles, respectively. Excitation energies for 
the light particles are higher, and the minimum of the 
quasi-particle dispersion shifts toward zero momentum. 



We also compare the average of the light and heavy 
particle dispersion relations to the dispersion obtained 
in the equal-mass case. The two are very similar, much 
as they would be in standard BCS theory. The super- 
fluid transition temperature would nevertheless decrease 
with increasing mass ratio, as the average of the mini- 
mum excitation energies in each branch is significantly 
lower than the minimum energy at equal masses. We 
find this average (calculated by subtracting p) to be 
T)(r = 6.5) = 0.38(4), in comparison to the calculated 
value r)(r = 1) = 0.50(5) in the equal-mass case. The in- 
dividual spectra are also important. For example, radio 
frequency response experiments, which have been used 
to explore the gap in the equal-mass case 14|, could be 
designed to be sensitive to the individual light and heavy 
particle dispersion relations. 

Away from equal populations, we explore the phase 
diagram at zero temperature by considering normal and 
gapless supcrfluid states using from 60 to 90 particles. 
For calculations of the normal state, the trial wave func- 
tions dictating the nodes are taken from free-particle 
Slater determinants with filled-shcll configurations in pe- 
riodic boundary conditions. In Fig. [2] we plot the ground- 
state energy versus the polarization P = (N/ l —Ni)/(Nh + 
Ni), in units of Ep£. 

The circles are QMC calculations of the normal state 
and the curve is a simple polynomial fit to the normal 
state results as a function of polarization. The poly- 
nomial fit is explicitly tied to the free-particle results at 
P = ±1, and to the binding energy B^m of a single heavy 
(light) particle in a sea of light (heavy) particles. With 
the majority particle number Ni^ and the simulation 
volume L 3 constant, we find B h = Q.36E l F = 0.99(5)E F lr 
and Bi = 2.3E F = 0.97(5) Eg r at r = 6.5. With this def- 
inition the equal mass binding B = O.GEp. At constant 
total densities, these results correspond to Bh = 0.76E l F 
and Bi = 2.7 E F , in rough agreement with results in Ref. 
[22| . By fitting the dispersion of single impurities, we 
find m*/mi = 1.3 and mt/mh = 1.0. 

We find an energy minimum near a polarization of 0.5, 
corresponding to a ratio of 3:1 heavy to light particles. 
The small value of the energy indicates a possible collapse 
of the normal state at a mass ratio smaller than that 
found in 3-body calculations, where Efimov states and a 
collapse begin at a mass ratio of 13.6 3, 24 1. At uni- 
tarity with this finite range potential we observe collapse 
(large negative energies and several particles within the 
interaction range) before r = 12. At a mass ratio r = 6.5, 
we find that the energy decreases quickly with interac- 
tion strength, reaching zero at hp a w 2.5. These few- 
particle correlations may increase loss rates and limit the 
effectiveness of standard cooling techniques which sweep 
from the Bose-Einstein condensation regime to the uni- 
tary regime. As the mass ratio increases, the minimum in 
energy will shift toward higher polarizations. It would be 
very interesting to examine this evolution and the associ- 
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FIG. 2: (color online) Equation of state of a heavy-light 
fcrmion mixture as a function of the polarization at constant 
density. Shown are results for the normal state (circles) as 
well as for the superfluid state (squares). The solid line is a 
polynomial best fit to the QMC results for the normal state. 



ated Fermi or Bose condensates of odd and even clusters 
of fermions. 

We also consider the possibility of polarized superflu- 
ids; a simple case is the gapless superfluid where unpaired 
particles are placed at the minimum of the dispersion 
curves in Fig. [T] The energies for the gapless superfluid 
state are shown as squares in Fig. [2J Over a range of po- 
larizations P < 0, we find the polarized superfluid has a 
significantly lower energy than the homogeneous normal 
state at the same density and polarization. 

The results displayed in Fig. [5] can be used to deter- 
mine the stability of these phases. We use a polynomial 
fit to the normal state to calculate the critical concentra- 
tions of heavy and light particles and possible first-order 
phase transitions that occur between the superfluid and 
the normal states at finite polarization. This is illus- 
trated in Fig. where we plot the energy normalized 
to the free-particle energy of the majority species at the 
relevant density to the 3/5 power as a function of the 
concentration x' = n^/ni for the majority light parti- 
cle case and x = ni/rih for the majority heavy particle 
case. The transition points can be found by equating the 
pressures and chemical potentials of the normal and su- 
perfluid states. For the case of majority light species, the 
equilibrium concentration of heavy particles is extremely 
small, x' c = 0.02(2), indicating equilibrium between a su- 
perfluid and a nearly fully polarized sea of light particles. 
For the majority heavy case the critical concentration is 
x c = 0.49(5), near the concentration found for the equal 
mass case: x c = 0.44 [l(|. The transitions indicated in 
the figure are calculated from the polynomial fits and 
indicated by dashed lines following the tangent construc- 
tion used in Ref . 11 1 . 
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FIG. 3: (color online) Equation of state as a function of the 
concentration x' (left panel) and x (right panel). Normal (cir- 
cles) and superfluid states (squares) are shown. The dashed 
lines are the coexistence lines between the normal and the 
unpolarized superfluid states. The resulting critical concen- 
trations are x' c = 0.02 and x c = 0.49. The dotted curves rep- 
resent non-interacting impurities with the calculated binding 
energies and effective masses. 



The polarized superfluid results are also shown in Fig. 
[3] Over a range of polarizations P < these states are 
very close to stability with respect to the phase separated 
normal state and unpolarized superfluid. It is possible 
that further generalizations of the trial states, for exam- 
ple by considering inhomogeneous polarized superfluids 
like LOFF states, would lower the energy and provide a 
stable polarized superfluid at zero temperature. 

The pressures and chemical potentials calculated for 
the superfluid and normal states can also be used to ex- 
amine what happens in a harmonic trap. Keeping the 
chemical potentials 1^1^ fixed and choosing the state of 
highest pressure with local chemical potentials fJ-h(i)( r ) = 
^h(i) ~~ Vh(i){.r) one can calculate, within the local-density 
approximation, the density for each species in the trap. 
In general, the trapping potentials of the two species are 
unequal; for this analysis we assume harmonic potentials 
with a strength raf t io\j1 for the heavy particles equal to 
twice that of the light potential strength miu>f /2, similar 
to that of a recent experiment on 6 Li - 40 K mixtures 17 1. 

In Fig. [l]we plot the local polarization as a function 
of scaled radius for various total polarizations P to t = 
(N h - Ni)/(N h + JVj) where N h and Nt are the total 
number of heavy and light particles in the trap. Curves 
are shown for P to t = -0.4, 0, 0.4, and 0.8. The radii 
are scaled in each case so that within the local density 
approximation the density falls to zero at r sc = 1. In 
this plot we assume that the polarized superfluid is un- 
stable at T = 0. For large total negative polarizations, 
the equilibrium configuration is an unpolarized (P = 0) 
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FIG. 4: (color online) The polarization of a trapped system 
as a function of the radius (scaled so that at r ac = 1 the 
density goes to 0). Shown are curves for four different values 
oiP t „t = (N h -N l )/(N h + N l ). 



superfluid in the center and a nearly fully polarized sea 
of light particles in the exterior. At a finite temperature 
the polarized superfluid state at P < would appear, 
similar to what happens in the equal mass case. Because 
of the small energy differences, we expect finite temper- 
ature effects to be even more important here. 

Near zero total polarization, in contrast, three distinct 
regions exist. In the center a normal state is favored with 
a polarization near P = 0.5; this is the lowest-energy 
normal state of the system described earlier. At larger 
radius, there is a shell of unpolarized superfluid, and then 
again in the exterior a region of nearly all light particles, 
or potentially a polarized superfluid. At this large mass 
ratio, this unusual configuration is actually lower in en- 
ergy than a homogeneous superfluid everywhere. Finally, 
for very large total polarizations the system is normal ev- 
erywhere, with the polarization smoothly increasing from 
the center of the trap as the radius increases. This is 
again analogous to what happens in the equal-mass case 
for a large total polarization. It would be interesting to 
confirm this new structure experimentally, to explore the 
stability and structure of the polarized superfluid, and to 
determine its evolution with mass ratio. 

In summary, we performed QMC studies of heavy-light 
fcrmion mixtures at unitarity. We find that the ground- 
state energy of the superfluid and the average quasi- 
particle dispersion agrees closely with the superfluid with 
equal masses and the same reduced mass m r , as predicted 
by BCS theory In contrast, the system at finite polariza- 
tion is very different from the equal-mass case, resulting 
in significantly different profiles of trapped systems, even 
for the case of equal numbers of heavy and light parti- 
cles. Polarized supcrfluids are very near to stability for 



more light than heavy particles. In the majority heavy 
case, we find an energy minimum at approximately a 3:1 
heavy to light ratio which evolves rapidly with mass ratio 
and interaction strength. The extra scale made available 
by the mass ratio produces a variety of fascinating new 
physical effects in cold Fermi atoms near unitarity. 
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Note added in proof— An analysis of the phase dia- 
gram of Fermi mixtures at unitarity was recently pub- 
lished [28| . It would be interesting to see this calculation 
repeated for the equation of state that we have obtained. 
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